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Abstract 

Cuprate superconductors have many different atoms per unit cell. A large fraction of cells (5- 
25%) must be modified ("doped") before the material superconducts. Thus it is not surprising 
that there is little consensus on the superconducting mechanism, despite almost 200,000 papers 
[TJ. Most astonishing is that for the simplest electrical property, the resistance, "despite sustained 
theoretical efforts over the past two decades, its origin and its relation to the superconducting mech- 
anism remain a profound, unsolved mystery [2]." Currently, model parameters used to fit normal 
state properties are experiment specific and vary arbitrarily from one doping to the other. Here, 
we provide a quantitative explanation for the temperature and doping dependence of the resistiv- 
ity, Hall effect, and magnetoresistance in one self-consistent model by showing that cuprates are 
intrinsically inhomogeneous with a percolating metallic region and insulating regions. Using simple 
counting of dopant-induced plaquettes, we show that the superconducting pairing and resistivity 
are due to phonons. 
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Since superconductivity requires coherent Cooper pairing of electrons, knowing what 
couples most strongly to electrons is absolutely necessary for understanding what causes 
high-temperature cuprate superconductivity. The resistivity, the Hall effect, and the mag- 
netoresistance are fundamentally measurements of the momentum dependent Fermi surface 
scattering rate, l/r(/c), that measures the strength of whatever is coupling to electrons. If 
the origin of the scattering rate and its temperature dependence are not understood, then 
most likely cuprate superconductivity is not understood either. 

The earliest resistivity (p) measurements [3] on cuprates found p to be approximately 
linear in temperature over the huge temperature range of 10 < T < 1000 K. A linear T 
resistivity is characteristic of electron-phonon (or electron-boson) scattering. Yet phonons 
are not believed to cause p, despite the fact that historically, the dominant scattering mech- 
anisms of electrons in metals have been phonons and impurities. There are two reasons 
for this conclusion: First, at low temperatures, the Bose-Einstein statistics of the phonons 
reduces the phase space for scattering, leading to p(T) ~ T 5 (the Bloch-Gruneisen law) [3]. 
The T 5 scaling should be observable for T < 0Dcb y c/lO where 0Deb ye is the Debye temper- 
ature (the characteristic energy of the highest energy phonons). Since 0Deb y e ~ 400 K in 
cuprates j5], phonon scattering is not compatible with the observed linearity. Magnons are 
also eliminated because G>Mag ~ 1500 K [6]. 

Second, the magnitude of p at high T for some dopings exceeds the Mott-Ioffe-Regel 
limit (MIR) [7j that occurs when the electron mean free path reaches the shortest Cu— Cu 
distance (~ 3.8 A). For cuprates, p(T) should saturate to Pmir ~ 1000 /if2-cm. Instead, p 
increases linearly right through pmir, leading to the conclusion that the normal state may 
not even be a typical metal (a non- Fermi liquid). 

Both of these conclusions are invalid when the crystal is intrinsically inhomogeneous, as 
we will show. Briefly, the first point implicitly assumes that phonon momentum is a good 
quantum number. The second underestimates Pmir by overestimating the density of charge 
carriers. For example, at optimal hole doping of x = 0.16 (we define x to be the number of 
holes per planar CUO2), we find the fraction 4 x 0.16 = 0.64 of the crystal is metallic leading 
to a pmir that is (1/0.64) = 1.56 times larger than the conventional estimate. p(T) remains 
below the larger Pmir up to the melting temperature. 

Recently, the doping and k- vector dependent scattering rate, l/r(x, k), has been extracted 
from a series of beautiful experiments by Hussey et al. [2J ISTfTO"] . In the first set of experiments 
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[U[9], the magnetoresistance (MR) of large single crystals of single layer Tl 2 Ba 2 Cu06+5 was 
measured as a function of the direction of an applied 45 Tesla magnetic field for T < 40K 
and large overdoping. They found that l/r(x,k) is the sum of three terms, l/r(x,k) = 
1/to(x) + Ai(x) cos 2 (2ip)T + A2(x)T 2 , where A\(x) and A 2 (x) are hole doping dependent 
constants and <p is the angle between the Cu— bond direction and the k- vector as shown 
in the top right corner of figure 6a. The first term, 1/to( constant that is sample 

dependent even at the same doping. The anisotropic linear T term is zero for fc-vectors along 
the diagonal (<p = 7r/4) and largest for k along the Cu— O bond directions (ip = 0, vr/2). 

In an elegant theoretical analysis, Hussey et al. [21 HU H2] showed that a scattering rate 
of the above form explained the resistivity, the Hall effect, and the MR for a large range 
of dopings, temperatures, and different cuprate materials. The authors suggested that the 
isotropic T 2 part of the scattering rate arose from electron-electron Coulomb scattering while 
the origin of the T term is unknown. Finally, the strong doping dependence of Ai(x) and 
A 2 (x) was noted, but not explained. 

In another experiment, Hussey et al. [21 [TO] used large magnetic fields to lower the T c 
(the superconducting transition temperature) of La2_a.Sr2.CuO4. They found that p(T) 
was the sum of doping dependent T and T 2 terms at low temperatures plus a sample 
and doping dependent constant. These results suggested that the 1/t(x, k) form found for 
Tl 2 Ba 2 Cu0 6+ 5 was applicable to all dopings and cuprate materials. 

The authors also found that the low temperature p evolved into linear T at high tem- 
peratures, as seen previously [13] ■ They concluded that the isotropic T 2 term was due to 
electron-electron Coulomb scattering and Mott-Ioffe-Regel saturation of the scattering rate 
at high temperatures was reducing this term down to linear T. 

The experiments described above may be summarized as follows [2j. The temperature 
derivative of the resistivity, dp(x, T)/dT, is of the form dp(x, T)/dT = oi(x, 0) + 2a 2 (x)T at 
low temperatures and dp(x,T)/dT = ai(x,oo) at high temperatures. ai(x,0) and ai(x,oo) 
are both doping dependent constants with ai(x, 0) < ai(x, 00) always, p smoothly interpo- 
lates between these two forms. In addition, «i(x,0) is approximately constant for x < 0.19 
and ai(„, 00) is approximately constant for x > 0.19. Surprisingly, the low-doping ai(x,0) 
constant and the high-doping ai(x, 00) constants appear to be equal (~ 1 pfl-cm/K) within 
the experimental error bars. A plot of the experimental data [2] for 0-1(0) and 0-1(00) for 
La 2 _~Sr~Cu0 4 is shown in figure 1 (for notational convenience, we suppress the x in Oi(0), 
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«i(oo), and l/r(k) from this point onward). 

Even for a fixed doping value, it is clear that the low and high T forms of l/r(k) put strong 
constraints on theory. An even stronger constraint is explaining the doping dependence of 
«i(0) and «i(oo) and their peculiar "crossover" at doping of x w 0.19, as can be seen in 
figure 1. 

In prior publications [HI [15], we showed that cuprates are comprised of a percolating 
metallic region and insulating regions. The metallic region is formed from 4-Cu-site square 
plaquettes in the CuC>2 planes that may overlap (share a Cu atom). These plaquettes 
are centered around the dopants (e.g., Sr in La 2 _ a; Sr a ,Cu04). The insulating regions are 
comprised of localized d 9 Cu spins that are coupled antiferromagnetically. 

By simply counting the volume and surface areas of these two regions as a function 
of doping, we explained the origin and doping dependence of the pseudogap (PG) [14J, 
the universal room-temperature thermopower, the STM incommensurability, the neutron 
resonance peak, and the generic cuprate phase diagram [T5]. Here, we show that cti(0) and 
ci!i(oo) are also explained by the same intrinsic inhomogeneity (figure 2). 

We will show that «i(0) is due to phonon scattering inside non-overlapping plaquettes and 
cti(oo) is due to scattering from all phonons. Therefore, ai(0) ex N^m and cti(oo) oc Nc u , 
where N^m equals the number of planar Cu sites inside non-overlapping plaquettes and Nq u 
is the total number of planar Cu sites. Dividing by the total number of charge carriers, 
N M , leads to ai(0) = C(N4 M /N M ) and ai(oo) = C(N Cu /Nm)- The best fit (see figure 1) is 
C = 0.904 /iO-cm/K. Figure la has one adjustable parameter, C. Figure lb is the ratio of 
ai(0)/ai(oo) and has zero adjustable parameters. The excellent fit in figure 1 using simple 
counting is the main result of this paper. 

Figure 2 shows a schematic of the distribution of plaquettes and how N±m , Nm, and Nq u 
are calculated. The insulating/antiferromagnetic (AF) regions (blue, green, and purple ar- 
rows) and the metallic region (yellow) are shown. The non-overlapping 4-Cu-site plaquettes 
are shown by black squares with the number "4" inside. The number of Cu sites inside these 
squares equals N±m- Nm is the number of Cu sites in the metal region (yellow) and Nc u is 
the total number of Cu sites (here, Nc u = 20 x 20 = 400). For x < 0.187, there is enough 
space to avoid plaquette overlap leading to N iM = N M , as seen in the x = 0.15 figure. The 
Methods section describes how the plaquettes are doped in figure 2. 

Prior to describing our detailed model for the anomalous normal state transport, we 
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summarize the cuprate doping phase diagram. All cuprate superconductors have CuC>2 
planes. When undoped, they are insulating Heisenberg spin- 1/2 antiferromagnets (AF) 
with a Neel temperature T N ~ 200 K and spin-spin coupling J dd m 0.13 eV [B]. The spin 
is localized on planar Cu atoms in a d 9 configuration with the single unpaired spin in the 
Cu d x 2_ y 2 orbital, where the x and y axes point along the Cu— O bonds. Doping of more 
than x ~ 0.05 holes per planar Cu destroys the bulk Neel AF phase and superconductivity 
occurs for T < T c where T c is the superconducting transition temperature. Above T c , the 
material is metallic. When x is increased past 0.27, superconductivity vanishes (T c = 0). 
The highest T c occurs for x ~ 0.16 [TH]. The band structure consists of approximately 2D 
Cu02 bands with Cu d x 2_ y 2 and O p a (p x and p y ) orbital characters that cross the Fermi 
level [T7j. The Fermi surface is hole-like and centered around k = (it /a, it /a) where a is the 
Cu02 unit cell distance (a w 3.8 A). 

In our Plaquette model [T4], doping creates holes in out-of-the-Cu0 2 -plane orbitals that 
are localized in square 4-Cu-site plaquettes centered at the dopant. This orbital character 
was shown by high quality ab-initio quantum mechanics calculations on La2_ a; Sr a; Cu04 su- 
percells [18J. These hole orbitals induce derealization of the 4 Cu d x 2_ y 2 and 8 O p a orbitals 
inside the doped plaquette. Each plaquette can be thought of as a tiny piece of metal (see 
Supplement). 

Figure 3 shows a Sr doped 4-Cu-site plaquette (as found in L^-zSr^CuO/t). In figures 3a 
and 3b the two degenerate out-of-plane states are shown with their respective Jahn- Teller 
atomic distortions (blue arrows for planar O atoms and black arrows for apical O atoms). 
In the cuprates, no long-range static distortion compatible with figure 3a or 3b is found. 
Thus the coupling between these out-of-plane hole states is strong, leading to a vibronic or 
dynamical Jahn- Teller state being formed instead. 

Figure 3c shows an instantaneous configuration of the out-of-plane holes in a single Cu02 
plane. The out-of-plane hole states are shown by the red "dumbbells". Coupling between 
adjacent plaquettes leads to correlation between the orientations of neighboring dumbbells. 
When plaquettes overlap as shown by the green squares in figure 3c, the orbital degeneracy 
of the out-of-plane holes shown in figures 3a and 3b is removed. Thus overlapping plaquettes 
have no dynamical Jahn- Teller state. 

In figure 3d, the correlation energy between neighboring dynamical Jahn- Teller states is 
shown for plaquettes inside a single Cu0 2 plane, e p i anar , and for plaquettes in adjacent Cu0 2 
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planes, e pe rp. In the Supplement, e perp is estimated to be ~ 3.6 x 10 5 eV ~ 0.42 K and 
e planar is on the order 0.01 to 0.1 eV (~ 100 - 1000 K). 

For temperatures greater than ~ 1 K, there is no correlation between dynamic Jahn- 
Teller states between planes. In analogy to the "dynamical detuning" proposed for electron 
transport normal to the CuC>2 planes [191 EQ], the lack of correlation of dynamical Jahn- 
Teller states normal to the plane disrupts the phase coherence of phonons between CuC>2 
layers. Hence, the phonon modes inside non-overlapping plaquettes become strictly 2D for 
temperatures greater than e pe rp. 

Since momentum is not a good quantum number for the 2D phonon modes inside the 
region of the non-overlapping plaquettes, these modes have the same character as those 
in amorphous metals. We show later that these amorphous 2D phonons lead to the low- 
temperature linear resistivity (the ai(0) term in figure 1). Also, ai(0) oc N^m because the 
number of amorphous 2D phonon modes is proportional to the total number of Cu sites 
inside the non-overlapping 4-Cu-site plaquettes, N iM . 

The phonon modes occurring "outside" the non-overlapping plaquettes are amorphous 
and 3D. The number of these modes is proportional to Nc u — Nam- These modes lead to 
the T 2 resistivity term, as we show later. 

A percolating pathway of plaquettes leads to a delocalized metallic Cu d x 2_ y 2 and O p a 
band inside the percolating swath (yellow region in figure 2). The undoped regions (not 
part of the doped metallic swath) remain localized d 9 Cu spins with AF coupling. Thus 
cuprates are intrinsically inhomogeneous with a percolating metallic region and insulating 
AF regions. 

Our proposed inhomogeneity appears prima facie to be at odds with the recent observation 
of quantum oscillations in heavily overdoped Tl-2201 cuprates [2TH23] . In the Supplement, 
we estimate the Dingle temperature arising from the inhomogeneity in our model and show 
that it is compatible with these measurements. 

The phonon modes that lead to superconductivity involve the planar displacement of 
the O atoms at the interface between the metal and AF regions (figures 4 and 5). These 
longitudinal-optical (LO) O phonon modes are the softened modes seen by neutron scattering 
[24H26] for the superconducting range of dopings for momenta along the (tt, 0) and (0, 7r) 
directions. They are different from the Jahn- Teller modes (figure 3a,b) that lead to the 
linear T term in p. Figure 6 shows how these modes lead to the observed d-wave (d x 2_ y 2) 
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superconducting gap. 

To complete the derivation of the results in figure 1, we must show that phonon modes 
inside the non-overlapping plaquettes produce the anisotropic cos 2 (2<^)T scattering rate and 
that the remaining phonons contribute the T 2 isotropic term. In a typical metal with full 
translational symmetry, the scattering rate for a metallic electron with momentum k is jl] 
1/rik) ~ J2k> ,q{fr/ Mu> q )q 2 (1 — cos 8)nB{oJ q )S(k' — k — q) where the sum is over k' on the Fermi 
surface and q is the phonon momentum. M is the nuclear mass, u q is the phonon energy, 9 
is the angle between k and k', and n^(u) = l/(e nW//fcsT — 1) is the Bose-Einstein distribution 
(ks is Boltzmann's constant). The delta function maintains momentum conservation. The 
(1 — cos#) term becomes T dependent at low T because restricts the scattering to 

small-angles {9 <C 1). 

All phonons can scatter k to k' in cuprates because of their intrinsic inhomogene- 
ity Hence, the (1 — cos 9) term in 1/r (A;) has no T dependence. In this case |27J, 
1/rik) ~ J2k' \{fr/Moj\)(k' — k) 2 n^(u\) where A is summed over all phonon modes and 
the delta function has disappeared. The number of excited phonon modes at temperature 
T < 0Dcb y e is ~ T d where d is the dimension of the phonon modes |28j. At high T, all phonon 
modes are excited. At low temperatures, 1/r ~ (l/u)T d ~ T d_1 , because hu ~ ksT . Our 
expression for 1/r is identical to the Debye- Waller factor [U [28] that appears in the resis- 
tivity of amorphous metals (see Supplement for details). At high T, 1/r ~ T because all 
phonon modes become thermally accessible and the occupation number for each mode, n^, 
is proportional to T. 

The 2D Jahn- Teller modes in figure 3 have d = 2, leading to 1/r ~ T. This linear 
scattering rate persists for T > e pcip ~ 0.4 K. The remaining 3D modes lead to 1/r ~ T 2 . 
Since the phonon modes are amorphous, l/r(k) is independent of k for the 3D T 2 modes. 
The 2D Jahn- Teller modes are derived from the vibronic distortions in figure 3. These modes 
raise the Cu orbital energy for two diagonal Cu atoms in a plaquette and lower the orbital 
energy for the other two Cu atoms. Hence, their scattering is modulated by an envelope 
function « cos 2 (2y?). Thus l/r(A;) ~ cos 2 (2ip)T. This modulation of l/r(k) does not change 
the total scattering around the Fermi surface. Instead, it merely redistributes its weight 
around the surface. Hence, the proportionality constant in front of the expressions for «i(0) 
and ai(oo) is the same. 

Thus the doping and temperature dependence of the resistivity, Hall effect, and mag- 
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netoresistance of cuprates is due to the intrinsic inhomogeneity of these materials. This 
inhomogeneity arises from the dopant atoms and leads to a percolating metallic region and 
insulating antiferromagnetic regions. The percolating metallic region is comprised of 4-Cu- 
site plaquettes that are centered around dopants. The plaquettes may overlap each other. 
We show that the ratio of non-overlapping plaquettes to the total number of Cu atoms equals 
the ratio of the low and high temperature linear T terms in the resistivity, ai(0)/ai(oo), 
with no adjustable parameters. We calculate this ratio and find very good agreement with 
experiment. The resistivity is shown to be due to phonon scattering. The special "crossover" 
doping, x ~ 0.19, for ai(0) and a^oo) occurs at the doping when plaquettes first begin to 
overlap. We calculate this value to be x = 0.187 with no adjustable parameters. The low 
doping value of cti(0) and the high doping value of cti(oo) are found to be exactly the same 
in our model. This is observed by experiments within the error bars. We find an excellent 
fit to the doping evolution of «i(0) and ai(oo) with one adjustable multiplicative constant. 
We also provide a physical picture for the origin of the observed softened LO O phonon 
mode and show that it leads to a d-wave superconducting gap. 

We conclude that the anomalous normal state transport in cuprates provides strong 
evidence for intrinsic inhomogeneity, metallic percolation, phonon superconducting pairing, 
and phonon normal state scattering. 

Methods 

We dope plaquettes using a simple model for the Coulomb repulsion of the dopants. Since 
the dopants reside in the metallic region, the repulsion is screened over a short distance. 
Thus we assume plaquettes are randomly doped without overlap. This doping is possible 
up to x = 0.187. For x > 0.187, it is impossible to add dopants without overlapping 
another plaquette. Above x = 0.187, the most energetically favorable location for a dopant 
is on a site that dopes three undoped d 9 Cu atoms and one previously doped Cu atom. At 
x = 0.226, it is no longer possible to locate three undoped d 9 Cu atoms that are part of 
a 4-site square. Hence, two undoped d 9 Cu's and two previously doped Cu's are doped up 
to x = 0.271. At this point, there are no adjacent pairs of d 9 spins. Single undoped d 9 Cu 
doping occurs up to x = 0.317. Above this doping, there are no remaining undoped d 9 Cu's. 
The green squares in figure 2 show all the plaquettes that overlap another plaquette (share 
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at least one Cu atom). The percolating metallic region has become the whole crystal. We 
have not computed dopings beyond this value (x = 0.317). For the results in figure 1, we 
doped a 10 3 x 10 3 square Cu0 2 lattice 10 3 different times. N iM /N Cu , and N M /N Cu values 
were obtained by averaging over the ensemble (see figures SI and S2). 
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FIG. 1: Doping evolution of the linear T coefficients of the resisitivity at high and low tem- 
peratures (o:i(CI) and Qi(oo)). (a) Experiment |2j versus theory. cci(0) = C(N&m/Nm) and 
cki(oo) = C(Ncu/Nm) where Nm is the total number of metallic Cu sites (yellow in figure 2), 
N^m is the number of Cu atoms inside the non-overlapping 4-Cu-site plaquettes (black squares 
with the number "4" in figure 2), and Nq u is the total number of planar Cu sites (see Figures SI 
and S2). The low and high temperature linear T scattering rates arise from phonons and are given 
by 1/t(0) ~ N^mT and 1/t(oo) ~ Nc u T, respectively. Dividing by the total number of charge car- 
riers, Nm, leads to ai(0) and ai(oo). For C = 0.904 ^O-cm/K, the RMS and maximum absolute 
errors are 0.017 and 0.17 /ifi-cm/K, respectively, (b), plot of the ratio ai(0)/ai(oo) = N^m/Nqu- 
C cancels out here, so there are zero adjustable parameters. For doping x < 0.187, N&m = Nm 
because there are no overlapping plaquettes. Hence, ai(0) = C. The start of plaquette overlap at 
x = 0.187 (vertical dashed line) leads to the sharp discontinuity in «i(0) there. Since Nm — > Nq u 
with increasing doping, ai(oo) at high doping tends to the same constant C, as observed, (a) and 
(b) are the main results of the paper. 
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FIG. 2: The distribution of plaquettes (non-overlapping and overlapping) on a 20 x 20 lattice in a 
single Cu02 plane. Each black, green, pink, blue, and red square represents a plaquette centered at 
a dopant (Sr in L^-xSr^CuO^). The corners of the squares are planar Cu sites. The O atoms are 
not shown. A Cu d x 2_ y 2 /O p a metallic band forms inside the percolating region of the plaquettes 
and is shown in yellow (see Supplement). The non-overlapping plaquettes (black squares with the 
number "4" inside) will have dynamic Jahn- Teller distortions as shown in figure 3. The number 
of these amorphous 2D phonon modes is proportional to the number of Cu sites inside these 
plaquettes, Nam- Nm is the total number of metallic Cu sites (yellow). Nc u is the total Cu sites. 
Here, Nq u = 20 x 20 = 400. The arrows are undoped d 9 Cu spins that are antiferromagnetically 
(AF) coupled. Plaquettes do not overlap for x < 0.187, leading to Nam = Nm (see x = 0.15 figure). 
For x > 0.187, overlap is unavoidable (green squares). See Methods section. The red squares are 
isolated plaquettes with no plaquette neighbors. These plaquettes lead to the pseudogap (PG) 
|14j . They are not connected to the percolating "metallic" swath. The 6x6 black square in the 
upper right of x = 0.19 is expanded in figure 3c to show the 2D Jahn- Teller phonon modes and in 
figure 4 to show the O atom distortions. The 4x4 square in the lower right is expanded in the 
Supplement to show the vicinity of an isolated plaquette (red square). 
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FIG. 3: Dynamic Jahn- Teller character of doped plaquettes and their atomic distortions, (a) and 
(b), the localized out-of-the-Cu02-plane hole character surrounding a Sr dopant in L^-xSr^CuO^ 
There are two degenerate states. The blue and black arrows show local distortions that stabilize 
these orbital states. A resonance of these two configurations (a vibronic or dynamic Jahn- Teller 
effect) occurs. These hole orbitals have Cu d z i and apical O p z character (vertical white-green 
orbitals) that induce the planar Cu d x 2_ y 2 and O p a orbitals inside this plaquette (shown in yellow 
in (d)) to delocalize over the 4-Cu-sites inside the plaquette. When the plaquettes percolate through 
the crystal (yellow region in (c) and figure 2), a planar Cu/O d x 2_ y 2/p a metallic band is formed 
inside the yellow region |14^ 18j (see Supplement), (c) Looking down at a 6 x 6 square in the CuO"2 
plane taken from the upper right corner of x = 0.19 doping in figure 2. The red "dumbbells" show 
the instantaneous orientation of the out-of-plane orbitals in the square. The degeneracy in (a) 
and (b) is split in the two overlapping plaquettes in the upper right (green squares). No dynamic 
Jahn- Teller states are formed there, (d) The coupling energy between neighboring "dumbbells" 
inside a Cu02 plane and between neighboring planes. The coupling energy between planes is very 
small, 

£perp ^ 3.6 x 10 ^ eV ro 0.42 K (see Supplement for the calculation). For T > ^perp ^ 

0.42 

K, there is no correlation of the dynamic Jahn- Teller states in adjacent Cu02 planes. Hence, 
the phonon modes inside non-overlapping plaquettes lose phase coherence between the planes and 
become strictly 2D I2U| . Within each Cu02 plane, these modes are amorphous (momentum 
is not a good quantum number). Amorphous 2D phonons lead to the low-temperature linear T 
coefficient, «i(0), as shown in the text. 
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FIG. 4: 3D picture of the distortion of the O atoms between the metallic (yellow) and AF (red) 
regions. These modes are shown in greater detail in figure 5 and lead to d-wave superconductivity 
as shown in figure 6. 
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FIG. 5: The AF/Metal interfacial O phonon modes. Fully expanded 6x6 black square in the 
upper right corner of x = 0.19 doping from figure 2. The interfacial O phonon modes (red arrows) 
are the origin of the observed softened phonon longitudinal optical (LO) oxygen modes seen by 
neutron scattering with wavevectors along the Cu-0 bond directions [(0, 0) to (tt, 0) and (0, tt)] for 
the superconducting range of dopings [231(26] . These modes lead to a d-wave superconducting gap 
as shown in figure 6. 



15 



Repulsive Everywhere „. 
No S-Wave Gap (n,n) 




:+U-V > 
Repulsive 



(71,0) 




Phase Cancels Coulomb 
D-Wave Gap Formed 









(7C,7t) 








i Repulsive 
V+(1/2)V>0 










-(1/2)V<0 
Attractive 


r 


(71,0) 












\ 












-k 2 | 



FIG. 6: Schematic diagram of the cuprate Fermi surface showing how a strong attractive phonon 
pairing along the Cu— O bond directions leads to the observed d-wave [d x 2_ y 2) superconducting 
gap. (a) and (b) show the hole-like Fermi surface centered at {ir/a, 7r/o) where a ~ 3.8 A is the unit 
cell distance (the T-point is occupied). The Cu— O unit cell in the lower left corner of (a) shows the 
orientation of the Brillouin zone relative to the crystal lattice. In both (a) and (b), a (k f, —k 1) 
Cooper pair in black is scattered by a phonon with momentum q\ or q<i to (k\ t, —k\ \) in green 
and (fo t> —^2 4) in blue, respectively. A Coulomb repulsion, (+U > 0), occurs for both q\ and q2- 
qi has an additional phonon- mediated attraction, V < 0, from the AF/metal interfacial O atoms 
in figures 4 and 5. For an isotropic s-wave gap, there is a net repulsion everywhere (upper left 
corner of (a)). Superconductivity cannot occur for this gap symmetry. For a d x 2_ y 2 gap (upper left 
of (b)), the average coupling around the Fermi surface of +U — (1/2) V integrates to zero leaving a 
net attractive — (1/2)V pairing along the Cu— O bond directions and a repulsive +(1/2V) coupling 
for q as (vr/a, ir/a). A d x 2_ y 2 gap can lower the energy. Thus the modes in figures 4 and 5 lead to 
the observed d-wave gap. The idea of a d-wave gap arising from the softened oxygen LO phonon 
modes along the Cu— O bond directions has been suggested previously |29j . In particular, Phillips 
|301 13"T] argued for softened modes occuring at the interface of the AF and metallic regions. What 
is new here, is that the physical origin of these phonon modes is provided along with its doping 
evolution. 
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Metal Plaquette Counting as a Function of Doping 

Figure SI shows the calculated number of Cu atoms inside the metallic region, Nm, and 
the number of metallic Cu atoms that reside in non-overlapping plaquettes, N^m- 111 figure 
2 of the main text, the number of Cu sites in the black squares with "4" inside is N±m and 
the yellow region of the same figure is Nm- The doping evolution of these two numbers 
determines the low and high T linear coefficients of the resistivity shown in figure 1 (ai(0) 
and ai(oo)). 

Figure S2 shows the doping evolution of the two different types of non-overlapping pla- 
quettes. Only the non-overlapping plaquettes inside the metallic regions (black squares with 
"4" inside in figure 2 of the main text) contribute to the low T coefficient of the resistivity 
(ai(0)). There are additional non-overlapping plaquettes in figure 2 outside of the metallic 
region. They are shown by red squares in the figure and are isolated plaquettes. The number 
and spatial distribution of these isolated plaquettes determines the pseudogap (PG) and its 
doping evolution as shown in reference [T4] . 

Cu/O d x 2 y 2/p cr Band Derealization Inside the Metallic Region 

Figure S3 shows the shift in the orbital energy of the metallic Cu d x 2_ y 2 orbital and 
the neighboring O p a orbital due to the dopant induced localized out-of-plane hole. Since 
the d x 2_ y 2 orbital energy is lowered more than p a , derealization can occur, leading to a 
Cu/O d x 2_ y 2/p a band inside the percolating metallic region. Further details are provided in 
reference [T4] and its supplementary information. 

Ab-initio DFT Calculations on La2 x Sr x Cu04 Supercells 

Pure density functionals such as local density (LDA) [SI, S2, 17] and gradient-corrected 
functionals (GGA, PBE, etc) [S3] obtain a metallic ground state for the undoped cuprates 
rather than an antiferromagnetic insulator because they underestimate the band gap due 



2 



to a derivative discontinuity of the energy with respect to the number of electrons [S4, 
S5]. In essence, the LDA and PBE functional include too much self-Coulomb repulsion. 
This repulsion leads to more delocalized electronic states. Removing this extra repulsion 
is necessary in order to obtain the correct localized antiferromagnetic spin states of the 
undoped cuprates. 

The self-repulsion problem with LDA has been known for a long time [S6]. Very soon 
after the failure of LDA to obtain the undoped insulating antiferromagnet for cuprates, 
several approaches were applied to correct the flaws in LDA for La 2 Cu0 4 . Using a self- 
interaction-corrected method (SIC-LDA) [S6], Svane [S7] achieved spin localization with an 
indirect band gap of 1.04 eV, and Temmerman, Szotek, and Winter [S8] found a band gap 
of 2.1 eV. Using an LDA + U method, Czyzyk and Sawatzky [S9] obtained 1.65 eV. In all of 
these calculations on the undoped cuprates, an increase in out-of-plane orbital character was 
noted in states just below the top of the valence band. Calculations with explicit dopants 
such as Sr in La 2 _ :E Sr. r Cu04 were not done. 

Our ab-initio calculations [S10, 18] were performed using the hybrid density functional, 
B3LYP. Due to its remarkable success on molecular systems, B3LYP has been the workhorse 
density functional for molecular chemistry computations for almost 20 years [Sll, S12]. For 
example [Sll], B3LYP has a mean absolute deviation (MAD) of 0.13 eV, LDA MAD = 3.94 
eV, and PBE MAD = 0.74 eV for the heats of formation, AH f , of the 148 molecules in 
the extended G2 set [S13, S14]. B3LYP has also been found to predict excellent band gaps 
for carbon nanotubes and binary and ternary semiconductors relevant to photovoltaics and 
thermoelectrics [S15, SI 6]. 

The essential difference between B3LYP and all pure density functionals is that 20% 
exact Hartree-Fock (HF) exchange is included. B3LYP is called a hybrid functional because 
it includes exact HF exchange. The HF exchange removes some of the self-Coulomb repulsion 
of an electron with itself found in pure DFT functionals. A modern viewpoint of the reason 
for the success of hybrid functionals is that inclusion of some exact Hartree-Fock exchange 
compensates the error for fractional charges that occur in LDA, PBE, and other pure density 
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functional [S17]. The downside to using hybrid functionals is they are computationally more 
expensive than pure density functionals. 

Our B3LYP calculations reproduced the experimental 2.0 eV band gap for undoped 
La2Cu04 and also had very good agreement for the antiferromagnetic spin-spin coupling, 
Jdd = 0.18 eV (experiment is ~ 0.13 eV) [S10]. We also found substantial out-of-plane apical 
O p z and Cu d z 2 character near the top of the valence band in agreement with LDA + U 
and SIC-LDA calculations. 

We also performed B3LYP calculations on La2_a.Sra.CuO4 for x = 0.125, 0.25, and 0.50 
with explicit Sr atoms using large supercells [18] . Regardless of the doping value, we always 
found that the Sr dopant induces a localized hole in an out-of-the-plane orbital that is 
delocalized over the four-site region surrounding the Sr as shown in figures 3 and S3. This 
result is in contrast to removing an electron from the planar Cu d x i^ y ijQ p a as predicted 
by LDA and PBE. 

Our calculations found that the apical O's in the doped CuOg octahedron are asymmetric 
anti-Jahn- Teller distorted. In particular, the O atom between the Cu and Sr is displaced 
0.24 A while the O atom between the Cu and La is displaced 0.10 A. XAFS measurements 
[SI 8] find the apical O displacement in the vicinity of a Sr to be pa 0.2 A. 

Compatibility of Our Proposed Inhomogeneity and the Quantum Oscillations in 
Overdoped Tl-2201 

De Haas-van Alphen quantum oscillations have recently been observed in heavily over- 
doped Tl-2201 [_D - [2"3"] . These authors concluded that there in no microscopic inhomogeneity 
in the Cu02 planes over a minimum of 1200 A. Hence, these experiments appear to contra- 
dict our model of microscopic inhomogeneity in the Cu02 planes. 

Here, we estimate the amount of scattering to be expected from the inhomogeneity in 
our model and show it is compatible with the data in Rourke et al [21J. Rourke et al 
find quantum oscillations in 3 of the approximately 100 samples tested. Two of samples 
showing oscillations had doping at « 0.31 doping and the remaing one had doping at ~ 0.26. 
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Significantly, no samples with lower doping showed quantum oscillations. 

Using the doping methodology described in the Methods section of our paper, we find 
that 4.5% of the planar Cu atoms are localized d 9 spins at 0.26 doping and 0.6% of the 
planar Cu atoms are localized d 9 spins at 0.31 doping. At these two dopings, only isolated 
d 9 spins remain in our model. Below « 0.26 doping, some of the d 9 Cu spins neighbor other 
d 9 Cu spins. These d 9 Cu atoms are no longer isolated. 

In his classic book, Shoenberg [S19] states that Dingle temperatures much above 5K 
smear out quantum oscillations so that they are not resolvable (page 62). Shoenberg also 
states (page 62 and Table 8.1) that every 1% impurity typically adds 10 — 100 K to the 
Dingle temperature. Thus it is possible the ~ 0.6% "impurities" we expect for the 0.31 
doped samples may still lead to observable quantum oscillations, but the 4.5% impurities at 
0.26 doping appears to be too large to be compatible with experiment. 

The flaw in the above simplified analysis is that a localized Cu d 9 "impurity" is not a 
typical impurity. The charge distribution around a metallic and localized d 9 Cu are virtually 
identical in the Cu02 planes in cuprates. For both Cu atoms, the Cu has approximately 
+2 charge with a hole in the d x 2_ y 2 orbital. The d 9 Cu impurity is not a charged impurity. 
Instead, it is a neutral impurity with exactly the same sized valence orbitals as the metallic 
Cu atom. Since the impurity scattering cross section (and hence the scattering length) scales 
as the square of the difference in the potential between these two types of Cu atoms, the 
scattering can be up to two orders of magnitude smaller than a typical impurity, as we show 
below. Thus 4.5% d 9 Cu impurities may lead to the same amount of scattering as 0.045% 
atomic substitutions of different atoms. At 10 — 100 K Dingle temperature per 1% impurity, 
this leads to a Dingle temperature of 0.45 — 4.5 K for the 0.26 doped sample. 

We can use the Friedel sum rule [I] to quantify this physical argument. The Friedel sum 
rule incorporates the response of the metallic electrons around an impurity to screen out the 
impurity charge distribution and maintain a constant Fermi energy throughout the crystal. 
The relevant Cu orbitals near the Fermi level are d x 2_ y 2, d z 2, and 4s with symmetries Bi g , 
Ai g , and Ai g , respectively. The charge difference between the impurity Cu atom and the 
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metallic Cu atom is Z = 0, leading to Z = = -[6a. + 6b] where 6b is the d x 2_ y 2 phase 
shift at the Fermi level and 6a is the combined d z 2 and 4s phase shifts at the Fermi level. 
Since there is predominantly d x 2_ y 2 character at the Fermi level, the d z % and 4s phase shifts 
should be small, 6a ~ 0.1. 

The Friedel sum rule therefore forces the d x 2_ y 2 phase shift to be small, 6b ~ 0.1 rather 
than a 6b ~ 1 for a "standard" impurity arising from atomic substitution. 

Since the scattering cross section scales as sin 2 6b [I], the scattering cross section may 
be two orders of magnitude smaller than a typical impurity. Thus we have shown that the 
scattering length expected from the inhomogeneity of our plaquette model leads to Dingle 
temperatures from 0.45 — 4.5 K at 0.26 doping and 0.006 — 0.6 K at 0.31 doping. 

Rourke et al measure the attenuation of the oscillations due to the Dingle temperature, 
but they do not quote their Dingle temperature. We can extract their Dingle temperature 
in two ways. 

First, the Dingle temperature, xjj, is given by 

k B x D = j~ (1) 

where ks is Boltzmann's constant and r is the scattering time. Multiplying the top and 
bottom of the right-hand-side of this equation by the Fermi velocity, Vp, we obtain 

7 hvp , . 

kBXD = 2ri Q (2) 
where 1$ is the scattering length. 1$ = 400 A in Rourke et al. Vp ~ 1 — 3 x 10 7 cm/s in 
cuprates [11]. For vp = 2 x 10 7 cm/s, xp> — 6.1 K. 
Second, we may use 

/ \ 1/2 

m e i ec \ /2/ie\ (2KF \ fl 



\m therm J \2tt / \ e J \l J 
where m t h eTm /m e i ec is the ratio of the band electron mass measured by the temperature 
dependence of the oscillation amplitude (~ 5 in Rourke et al), fis is the Bohr magneton, F 
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is the oscillation frequency (1.8 x 10 4 Tesla in Rourke et al). This leads to xd = 5.2 K. 

Therefore, estimates of the scattering from our plaquette induced inhomogeneity lead to 
a Dingle temperature compatible with observations. 

In our model, below 0.26 doping, there are d 9 Cu atoms that are no longer isolated. 
These Cu atoms form antiferromagnetic regions in the material. We anticipate that the 
scattering from these "non-isolated" regions will lead to greater scattering than isolated d 9 
Cu atoms. This may explain why no quantum oscillations were found for dopings less than 
« 0.26. 

Despite the calculations done here, further work needs to be done to obtain precise 
numbers for the scattering cross section from an isolated d 9 Cu spin. What we have shown 
here is that it is most likely that the inhomogeneity we propose in this paper does not 
contradict the quantum oscillations found in heavily overdoped Tl-2201. 

The Debye- Waller Factor and Amorphous/Glassy Metals 

For liquid metals, Ziman and Faber [S20, S21] showed that the T dependence of the 
structure factor, S(K), evaluated at K = 2k f, where fcp is the Fermi momentum determines 
the T dependence of the resistivity. Their resistivity expression is given by [S20, 28], 

Here vp is the Fermi velocity, t(K) is the scattering matrix element, and Q is the volume. 

The integral for p is dominated by the maximum momentum transfer across the Fermi 
surface, or 2kp. Thus the T dependence of the structure factor, S(2kF,T), determines the 
temperature dependence of p. 

In amorphous/glassy metals, the structure factor is dominated by the "zero-phonon" 
elastic term that is given by S(K,T) = S(K)e~ 2W W where S(K) is the static structure 
factor and W(K) is the Debye- Waller factor arising from the nuclear motion at temperature 
T [S22, S23, S24, S25]. Here, K = 2k F . Using S(K,T) « S(K)[1 - 2W(K)\, we conclude 
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that the temperature dependence of the W(K) determines the temperature dependence of 

Pin 

The Debye- Waller factor is given by W(K) = (1/2) Y^ q \K ' U q \ 2 where the sum is over 
all phonon modes q and U q is the amplitude for the q mode. U q can be determined from 
(l/2)Mojq\U q \ 2 ~ {l/2)n(uj)Tiuj q pi] where M is the nuclear mass, hu q is the energy of 
the phonon mode q, and n q is the Bose-Einstein occupancy of the phonon mode, n q = 
IfteP*"* - 1). From this expression, we find the form of the the Debye- Waller factor to be, 

W(K) ~K 2 [ ^d 3 q. (5) 

J UJ q 

The expression for W(K) has exactly the same form that we obtained from the matrix 
element argument in the main text |27j . 

The Bose-Einstein occupation, n q , cuts off the integral at ~ ksT ', where ks is Boltzmann's 
constant. Since u q = cq, where c is the speed of sound, W(K) ~ T 2 . For 2D phonon modes, 
the integral is over d 2 q instead of d 3 q leading to W(K) ~ T. 

For 2D, the Debye- Waller factor diverges as / du/u ~ logT as T — > 0. In fact, the 2D 
phonon modes are not strictly 2D for T < e pcrp ~ 0.4 K. Less than 0.1 K, these modes are 
3D leading to a finite value for the Debye- Waller factor. For T > 0.1 K, the Debye- Waller 
factor is linear in T leading to a linear resistivity. 

Estimating the Coupling Energy of Dynamical Jahn- Teller Plaquettes in Neighbor- 
ing CuO-2 Planes 

The coupling energy, e pe rp, between two dynamical Jahn- Teller plaquettes in adjacent 
planes is given by the Coulomb repulsion energy difference between the different orientations 
of the plaquette hole orbitals. The details of the energy difference expression are shown in 
Figure S4. 

The Thomas- Fermi screening length, A, is estimated from the expression [1], 

4 = 47re 2 iV(0) (6) 
A 
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where e is the electron charge and N(0) is the density of states. There is approximately one 
metallic electron per unit cell leading to a density of states of approximately one state per 
eV per unit cell. The volume of the unit cell is, Q ce u = a 2 d, where a = 3.8 A and d = 6.0 A. 
Therefore, N(0) « 87.0 eV^A" 3 and A = 0.69 A. 

Substituting these numbers into the expression shown in figure S4 leads to e pe rp = 3.6 x 
1(T 5 eV = 0.42 K. 

Oxygen Modes Around an Isolated (Red) Plaquette 

The red arrows in figure S5 show the softened phonon mode for the O atoms between the 
AF spins and the delocalized electrons inside the isolated plaquette. 
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FIG. 1: Supplementary Material (FIG. SI): Calculated Nm and N^m as a fraction of the 
total number of Cu sites, Nc u - The green curve is exactly the same as ol\ (0)/ax (00) in figure 
1. The error bars in this figure are smaller than the width of the lines. The algorithm used to 
compute these curves is described in the Methods section. 
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4-Cu-Site Metal and Isolated Fractions 
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FIG. 2: Supplementary Material (FIG. S2): Doping evolution of the metallic and isolated 
non-overlapping plaquettes. The green line is the number of Cu sites contained in non-overlapping 
plaquettes inside the metallic region per planar Cu, N±m /Ncui an d the red line is the number of 
Cu sites inside the isolated plaquettes per planar Cu, N4i so i a t ec i/Nc u , appearing as red squares in 
figure 2. The green line is equal to ai (0)/ai (oo) in figure 1 and is also plotted in SI. The error 
bars in this figure are smaller than the width of the lines. The isolated plaquettes are responsible 
for the pseudogap [2] and are not relevant for the resistivity. The black line is the total number 
of Cu sites in non-overlapping plaquettes, N^Totai = N^m + ^iisoiated- For x < 0.187, there is no 
plaquette overlap leading to N^Totai/Ncu = 4x. For x > 0.187, Nii so i ate d — > very rapidly, leading 
to N 4M -> N iTotal . 
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FIG. 3: Supplementary Material (FIG. S3): Delocalization of Cu d x i_ y ijO p CT electrons 
inside 4-Cu-site doped plaquette [HI [18]. The additional hole induced by substituting Sr for La 
does not go into the planar Cu d x i_yijO p a orbitals as is usually assumed. Instead, it is comprised 
of predominantly apical O p z and Cu d z 2 as seen in figure 3 of the main text. The planar Cu 
sites surrounding the Sr dopant are called doped sites. Planar Cu sites that are not in the vicinity 
of a dopant (only La nearby) are called undoped sites. The energy level diagram shows how an 
out-of-plane hole above the planar Cu sites leads to delocalization of the d 9 Cu spins and the 
formation of delocalized Cu d x 2_ y 2/0 p a states inside the plaquette. The left figure shows the 
energy levels that occur at undoped Cu sites where there are no holes in out-of-the-plane orbitals. 
In this case, the difference in energy between the doubly occupied Cu site and the planar O orbital 
energy, e x 2_ y 2 + U — e p , is large, leading to localization of spins on the Cu. In the right figure, the 
Cu orbital energy is reduced due to the missing electron in the apical O sites directly above the Cu 
atoms leading to the neighboring doubly-occupied O p a electrons delocalizing onto the Cu sites. 
When two plaquettes are neighbors, the delocalization occurs over all eight Cu sites. When the 
doping is large enough that the plaquettes percolate in 3D through the crystal, a "metallic" band 
is formed in the percolating swath and current can flow from one end of the crystal to the other. 
This metallic band carries the current in the normal state and becomes superconducting below T c . 
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FIG. 4: Supplementary Material (FIG. S4): Estimating the coupling between dynamic 
Jahn- Teller modes in neighboring Cu02 planes. The coupling, e pevp , is given by the difference in 
Coulomb repulsion energy between the blue and red "dumbbells", e(B, R), and the orange and red 
"dumbbells", e(0,R). The 1/R Coulomb repulsion is attenuated by the Thomas- Fermi metallic 
screening length (dervied in the text to be, A ~ 0.69 A). Since the charge in a dumbbell is divided 
equally between the two ends, a charge of (l/2)e is used for each interaction term (e is the electron 
charge). The Cu-Cu lattice spacing in the plane is 3.8 A and the distance between planes is 6.0 A. 
We find e DerD = 3.6 x 10~ 5 eV = 0.42 K. 
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FIG. 5: Supplementary Material (FIG. S5): Oxygen phonon mode around an isolated 
plaquette. The figure zooms into the 4x4 dotted black square in the lower right corner of the 
x = 0.19 doping in figure 2. The yellow here denotes the isolated plaquette (shown by a red square 
in figure 2). Yellow is used here to signify that there is derealization of the Cu/O d x 2_ y 2p a orbitals 
inside the plaquette. Since this plaquette is not connected to any other plaquette, these electronic 
states remain localized inside the plaquette. An orbital degeneracy at the Fermi energy of these 
electronic states leads to the pseudogap [J5] . The O atoms that reside between the localized spins 
in the AF region and the delocalized electronic states inside the plaquettes will distort and the 
energy of this mode will change. We show the phonon mode by red arrows. 
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